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We calculate the shot noise power in ballistic graphene using the kinetic equation approach based 
on the Keldysh technique. We find that the local energy distribution function obeys Poisson's 
equation, indicating a mapping into a diffusive metal system. We derive the conductance and noise 
including the long range Coulomb interaction to first order. We find that the shot noise increases 
due to interaction, leading to a frequency dependence. Furthermore, we find that the Fano factor 
at degeneracy is 1/3, the same as without the Coulomb interaction. 

PACS numbers: 73.50.Td, 73.23.Ad, 73.63.-b 
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is also known in diffusive metals. 
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Graphene, a two-dimensional honeycomb lattice of car- 
bon, is of considerable recent interest. It is now experi- 
mentally established^ that a great deal of properties of 
graphene can be understood in terms of noninteracting 
linear dispersing Dirac quasiparticles. A particular atten- 
tion focuses on ballistic transport. To identify the elusive 
long range Coulomb interaction we study here shot noise, 
which is a fundamental nonequilibrium quantity. 

Conductivity and shot noise S n have been experimen- 
tally measured in graphene^—. The conductivity for a 
wide graphene ribbon (the width of the ribbon exceeds 
the its length: W >> L) was predicted^— to exhibit a 
minimum a m i n at the Dirac degeneracy point, while the 
Fano factor F = S n /(2ea m i n V) had a maximum close to 
1/3 (V is the voltage). The conductivity near the neu- 
trality point is attributed to the fact that the current is 
mediated by evanescent rather than propagating modes. 
The transport by evanescent channel is suppressed by a 
factor exp(-kyL) (here k y is the momentum in the width 
W direction). However, in the case W » L the expo- 
nent is actually small (k y L << 1) for many channels, 
resulting in an universal minimal value for conductivity 
at the Dirac point. 

The minimum conductivity at the graphene Dirac 
point can be reached in clean disorder-free graphene. To 
achieve a such regime may be a problem even for ultr- 
aclean high mobility suspended graphene^ with a possi- 
bility of inhomogenous charge landscapeiSdue to charged 
impurities. The measurements^ result in a conductance 
near neutrality point larger than theoretical values, while 
other data£ do obtain a m i n . The origin of this discrep- 
ancy is unclear at present. 

The Fano factor has been analyzed theoretically in 
a number of work o 8 ' 11 ! 12 relying on the Landauer ap- 
proach using a distribution of transmission eigenvalues; 
this approach is limited to non-interacting systems. At 
low temperatures a universal minimum conductivity is 
found cr m in = (4/ir)e 2 /h as well as a universal Fano fac- 
tor F = S/2eI = 1/3. Interestingly, the F = 1/3 result 



The Landauer approach is difficult to generalize when 
electron-electron interactions are present. Here we ap- 
ply the kinetic equation method, which is based on the 
Keldysh technique. This method is an alternative that 
permits to handle the Coulomb interaction near the Dirac 
point. A second advantage of our method is that the en- 
ergy distribution function can be readily identified. We 
address zero-frequency shot noise; however, at low u> a 
frequency dependence of the shot noise will be included 
through the Coulomb coupling rcnormalization. 

Our results can be tested by two specific experiments: 
first is the probe of local energy distribution as function 
of the distance R from one of the electrodes, as done 
e.g. in disordered wires^. We predict that the distri- 
bution has two steps and it interpolates linearly with R 
between the Fermi distributions with voltages 0, V of the 
two electrodes, respectively. A second experiment tests 
the presence of the Coulomb coupling g by looking at 
the conductivity or at the shot noise at finite frequency 
uj. The correction to the noise varies [Eq. below] as 
g r ~ l/ln(<j) where g r is the renormalize d 14 ' 15 interac- 
tion at finite oj. 

The paper is organized as follow. In Sec. II we intro- 
duce the Keldysh action of the system and present the 
formula of noise as a variation in the partition function 
with respect to the quantum components of the vector 
potential. In Sec. Ill we consider the noninteracting bal- 
listic graphene near the Dirac neutrality point. For this 
purpose we develop a new approach based on a quasi- 
classical approximation: by gradient expansion we obtain 
the kinetic equation for the non-equilibrium distribution 
function. This kinetic equation is a diffusion type similar 
to the one which describes one dimensional dirty wires. 
In Sec. IV we consider the impact of long range Coulomb 
interaction and present the results for conductance and 
shot noise power. The main calculations which are re- 
lated to this section are given by Appendices (A and B). 
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II. 



KELDYSH ACTION 



We use the standard form for the Hamiltonian that 
describes graphene, being equivalent to QED in 2+1 
dimension^,™ 



ic 



Mij x hv F (d x + j^A x ) + i^hvpdy}^ (l) 
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Here i\) a are 4 component spinors corresponding to two 
incquivalcnt Dirac points and to two atoms in the unit 
cell, cr is a spin index and A x is the x component 
of the vector potential. The Dirac matrices satisfy 
the standard algebra {7^,7,,} = 2<? MI , with the explicit 
representation^ 



7o 



where a are Pauli matrices. 

We write the corresponding action in the Keldysh ro- 
tated form. We consider the multichannel limit when 
the width W of the graphene layer is bigger than its 
length L, i.e. the aspect ratio is W/L ^> 1. This al- 
lows us to replace the summation over channels by in- 
tegration over momentum. The graphene layer occupies 
the space < x < L, and we consider the fluctuations 
of the total current (integrated over W), therefore only 
the x-component of the electromagnetic field is needed. 
Thus, the total action as a function of the vector po- 
tential (its classic and quantum components) acquires a 
form S (A) = S (A) + S int 



So(A) 



J2 J dPxdtyaid- 1 - ^1xA x )^ (2) 



where A x = A x tq + t x A x ; A x and A q x are the classi- 
cal and quantum components of the vector potential, 
respectively^. Here t x is a Pauli matrix and To is the 
unit matrix, both act in the Keldysh space, and ip a be- 
comes an 8-component spinor including the Keldysh in- 
dices. The Green's function has a form 



GR s~<K 
a 

Gt 



(3) 



with each entry as a matrix in Dirac space. The classical 
field A x is included in G -1 , hence 

(G^)~ 1 (A) = ifhd t + ij x hv F (d x + ^-A x ) + i^hvpdy 

(4) 

The unscreened long range Coulomb interacting is 
given by the part of the action which in the rotated 
Keldysh basic acquires a form (we set vp = fi = 1) 
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§X)/ d 2 xd 2 x'dt(^(xt)- / Q T i; ry ( X t)) i — 
(^ ff /(x't)7°r ai ^'(x'i)) (5) 



To first order in the interaction the partition function 
Z(A) = J D(ipip)exp(iS(A)) becomes Z(A) = Z (A)[1 + 
i < Si„t > (A)}, where Z (A) corresponds to the action 
S (A). 

The current-current correlation function can be ob- 
tained by taking the second derivative of the source- 
dependent partition function with respect to the quan- 
tum component of the vector potential: 



S n (t,t') = 



dR 
2L 



d 2 r 



S 2 InZ(A) 



SA%(xt)SA%(xV) 



(6) 



here R = (x + x')/2 and r = x x'. The leads serve 
as reservoirs of equilibrium electrons, and we take two 
arbitrary sections x, x' of the graphene area. A similar 
approach was undertaken for calculation of the shot noise 
in dirty wires^. 



III. NONINTERACTING BALLISTIC 
GRAPHENE 



In this case the noise power acquires the form 

S n0 (t,t') = «y^y ^ j drTr[G a {*t*!t') 

j x T x G a {x.'t'x.t)j x T x ] (7) 

where a a = W/L is the aspect ratio. In Eq.([7]) we 
perform the trace in the Keldysh space 

S n o(t,t') = ajj^f^rj drtr[G R a {*tx!t') lx G*{x!t'*t) 

(8) 



+G^(xtxY)7*G^(xYxi) 
+Gf(xfecY)7xG CT : (xYxi)]7 :c 



This involves Keldysh Green's function which has the 
standard parametrization 



G = G H F - FG 



(9) 



and satisfies Dyson's equation. The matrix function F 
is the nonequilibrium distribution function^. In equi- 
librium (V = 0) the Fourier-transform of this function is 
/o(e) = tanh(e/2T) (for fcrmions) . For graphene with the 
Dirac hamiltonian, F is block diagonal matrix with the 
two blocks F- and F+, each a 2x2 matrix. However, the 
diagonal elements are dominant at the Dirac point since 
the off diagonal elements are proportional to the energy 
deviation from the degeneracy point, i.e. vanish at the 
Dirac point. For the ballistic transport regime Dyson's 
equation is reduced to an equation for the matrix func- 
tion F which further can be transformed by performing 
gradient expansion into a kinetic equation. In the clean 
limit we have collision-ness limit for this kinetic equation 



F(G A )~\A) - (G R )~\A)F = 



(10) 
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The vector potential A x can be eliminated from this 
equation by a gauge transformation 

F(xt,xV) = C/(xt)F(xt,x't')[/ t (x't') (11) 

where U(xt) = exp[ze f Q dx'A x (x't)] and the block en- 
tries for F arc then F_ and F + . However, the boundary 
conditions for F(xt,x't') will be modified to include the 
phase factor—. By a standard procedure 1 ^ we obtain the 
first kinetic matrix equation 

[(_ ^ + t)F - + ^2-^+P^ F - a - + Py F - a y^ - 



cr x {- 







2dR 

Px)F + +p y a v F + ] = 



(12) 



Here r = (t + t')/2 and Fourier transform on the dif- 
ferences x — x', t — t' has introduced energy and mo- 
mentum variables e,p x ,p y . The second equation follows 
from Eq. (fl"2|) by replacing — > F+,- and <T x ,<r y — > 
—a x , —cr y . 

We can solve the kinetic equations independently for 
the F± block functions. It is easy to recover the rela- 
tion -F-(e) = F+(— e). Remarkably, we obtain a simple 
equation for diagonal function F d = (Fl 1 + F 22 )/2 which 
is, as we mentioned, of principal importance for noise 
calculation near degenerate Dirac point (e — > 0) 



d 2 F d 
dR 2 



k 2 F d = 



(13) 



where k = 2e^Jl + p 2 /(p 2 — e 2 ). Recalling that the 

boundary conditions for F include the phase factor we 
have at the boundaries 

F d (R = 0,tt') = /o(t-t') 
F d (R = L,t,t') = / (t-t')exp(-^(t)+^(O) 

(14) 

where the Fourier-transform of fo(t — t') is /o(e) and the 

phase <p(t) = Jq dx'A x (x't). For the constant applied 
bias V we can write F d (R, t, t') = F d (R, t — t 1 ). Focusing 

on the Dirac point, Ea.([i"3j) is reduced to = 0. The 
solution of this equation is 



F d (R,t,t') R R 



or in Fourier transform 
F d (R,e) 



(l-f)/o(e) + j/o(e 



(15) 



eV) (16) 



This two step distribution can be tested experimentally 
at intermediate positions < R < L 

In the next step we apply the transformation Eq. (|lll) ] 
to the each Green's function (GF) in the formula for the 



noise power ([7]), use representation (0 and notice that 
the energy integration over the product of GF of types 
G A G A or G R G R gives a zero result. The traces of re- 
tarded and advanced GF in the formula for noise are 
independent of R. Therefore, we can perform a direct 
integration in Eq.Q- The zero frequency noise becomes 

Sno = « J deK(e)[l-f 2 (e) + ^(2f 2 (e)- 

-f 2 (e-eV)-f (e)f (e-eV))] (17) 
where we define 

K ^ = ^TT,j -0ftr(lmG?(pe)j x ImG R (pe) 7x ). 

(18) 

The retarded GF has a form 



G*(pe) = 



(e + z£)7q - 7 ■ p 
(e + iS) 2 — p 2 



hence near the Dirac point ImG^'(pe) « — joS(p 2 + S 2 )" 1 
and we get (in standard units) n(e) = <J m in- 

Let us consider two limits of Eq. (fTTj) . the first is the 
equilibrium noise T > eV — > 0, 

S n o = " J de«(e)[l - /g(e)] = 4Taa min . (19) 

The shot noise corresponds to the limit V > T — > 0: 
a 



S n0 - 3 I de K (e)f (e)(f (e)-f (e-eV)) (20) 

Completing the integration we get S n o = 2e\V\aa m i n /3, 
hence the Fano factor is F — 1/3. 



IV. COULOMB INTERACTION 

The first order contribution to the noise due to 
Coulomb interaction ([5]) is presented by two topologically 
different sets of diagrams: one (Fig. la) is the self-energy 
contribution, the other (Fig. lb) is the vertex type dia- 
gram. We write them in terms of matrix GF's (in Dirac 
and Kcldysh space) . The detail calculations of these di- 
agram are given in Appendix A. Here we presents the 
result. In the case of equilibrium noise Coulomb contri- 
bution is 

S g = 0.76aTa min g r 

We note that this value of noise comes out of the nondi- 
vergent part of vertex diagram (Fig. lb) while logarithmic 
divergences in both self-energy and vertex diagrams can- 
cel. 

Here we introduce the Coulomb renormalized coupling 
g r that defines the interaction contribution to the noise. 
It is known to decrease as a function of frequency and flow 
to zero 14 ' 15 in the limit ui — >• 0, i.e., g r = 87r/ln(A/w), 
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71). Thus in standard units we have 



2e|VI 



5 g = a — - — 0.l9a min g r 



FIG. 1: A set of diagrams describing the first order correc- 
tions to the noise. The dotted lines represents the current 
vertex, the solid lines stand for retarded or advanced Green's 
functions, the solid lines with a triangle represent Keldysh 
functions, and the dashed line refers to the Coulomb interac- 
tion. An additional set of diagrams involves only one Keldysh 
and three retarded or advanced Green's functions. 



A 




V 

FIG. 2: The two loops diagram which describes the first order 
corrections to the noise. The lines are the same as in Fig.l. 
The additional diagrams of this kind are those which replace 
the loop with two Keldysh Green's function (like left one in 
Fig. 2) by retarded and advanced Green's functions. 



where A is an ultraviolet cutoff. This uj dependence 
can be probed experimentally to identify the effect of 
Coulomb interactions. 

For finite frequency w, however, there are additional 
diagrams ( like one showed in Fig. 2 ), which connect two 
electron loops by an interaction line, in analogy with the 
dirty wire case^S. In Appendix B we show that the shot 
noise contribution of this diagram for low frequency and 
uj > q is proportional the the small transmitted momen- 
tum q, i.e. vanishes as q — > 0. Therefore, we can neglect 
the diagrams with two electron loops. 

The renormalized g r is finite at w = if weak disorder 
generated by ripples is present. An attractive line of fix 
points appears then in the flow diagra m 14 ' 15 leading to 
a finite term in the interacting part of the noise even at 
uj = 0; however, other disorder corrections to the noise 
may arise, which are not addressed here. 

Next we consider the implication of Coulomb interac- 
tion in graphene on the zero frequency shot noise, i.e., 
V > T. We take the limit T — > and study the case 
of small voltages (we keep only terms linear in V). We 
substitute the representation © for Keldysh GF and the 
solution (|15p into Eqs. (| A3|) and (|A4|) . We find a common 
factor for both the self energy and for the vertex part of 
the noise, which at T —> coincides with the term in the 
square brackets of Eq. (fT7f (we note that an integration 
which involves only one Keldysh Green's function results 
in Fd — > 1 at T — > 0). This factor replaces the value 
1 — /q in the expressions for I2 (p) (|A6|) and in J\ (q) (Eq. 



(21) 

To find the Fano factor we extract from the equilib- 
rium noise the conductivity and, by this calculate the 
Poissonian noise: Sp = 2el = 2e\V\aa m i n (l + 0.19g r ). 
Comparing with the expression for the shot noise 

2e\V\ 

S = S n0 +S g = a—^-a min {l + 0.19g r ) (22) 

we find F = S/Sp = 1/3. Thus the unscreened Coulomb 
interaction precisely cancels and the Fano factor is inde- 
pendent of g. 

The value of Fano factor and its independence on in- 
teraction (at the first order) is in striking resemblance 
with analogous behavior in dirty wires 2 ^ when the limit 
of large inelastic length is considered. This similarity 
between clean and dirty systems is due to evanescent 
modes which define the current-current correlations in 
pure graphene at the Dirac point. The diffusive type 
equation for distribution function is common to graphene 
and to dirty wires. This distribution function preserves 
the Fano factor to the first order in g r . We remark that 
in other cases F does depend on the interaction to 1st or- 
der, e.g. in the interacting spinless resonance model and 
in the Kondo problem 2 ^. The Fano factor is also changed 
by strong interaction (in the so called hot electron limit) 
in dirty wires with a short relaxation length. However, 
this case is described by a quasiequilibrium distribution 
function and is beyond the applicability of perturbation 
theory^. In the case of ballistic graphene the hot elec- 
tron limit is unlikely since the interaction flows to small 
value at small uj. 

A few technical remarks are in order, related to the 
numerical coefficient of the coupling g. The exact value 
of this coefficient depends on the order in which uj and 
the finite inverse lifetime S are taken to zero. We take 
the dc limit uj — > while keeping 6 > (see discussion in 
Refsi 11 ' 12 ). The factor (1— f§ (e)) in the relevant integrals 
at T — > supports such a choice. For the self-energy di- 
agram we chose hard cutoff which restricts the fermion 
energy to the Dirac point 1 ^. In works that use the Kubo 
formula for the conductivity an additional rcgularization 
procedure for the self energy (Fig. la) is necessary, a pro- 
cedure that is a the subject of debat o 15 ' 23 " — . These dif- 
ferences may cause our numerical coefficient 0.19 to de- 
viate from that of Ref. [24| while being closer, though 
bigger, than the result 0.081 of Ref. [l||. 

V. CONCLUSION 

We apply the standard Keldysh technique and kinetic 
equation approach to calculate the shot noise power of 
clean graphene at the Dirac point for large aspect ra- 
tio W >> L. Considering e.g. a sample^ with W/L = 
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24, L = 200nm, the limitation of being at the vicinity of 
the Dirac point implies that V,ui,T < v F /L w 20meV. 
We suggest that tunneling experiments, as done in dirty 
wiresi^, can be done in this voltage range and can test 



the linear interpolation of Eq. (|16[) . i.e. a position de- 
pendence of a two-step energy distribution function. We 
have also found an interaction correction for the shot 
noise, varying as g ~ l/ln(w) (Eq. (|2"2"jl ). This coulomb 
effect may be detected by the frequency dependence of 
either the conductivity or the shot noise. We note that 
Vp/L is the ballistic flight time across the system, hence 
at cj <C Vp/L the noninteracting term is expected to be ui 
independent^. Therefore an observed u dependence in 
this range can identify the elusive Coulomb interaction. 
Since the same correction appears in the conductance 
and in the shot noise the Fano factor remains robust to 
interactions and persists being equal to 1/3. 
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Appendix A 

In this appendix we evaluate noise diagrams in Fig. 1. 
The self-energy diagrams can be written explicitely as 



e 2 iga 



2L 

G CT (xVyti) 



^2 / Tr [G <T {2tixt)j x T x G <J (xtx't')j x T x 

1 



l*-y| 

loT x G a {yt 1 zti)^ )] 



(Al) 



e 2 iga 
2L 

1 



St 



e 2 iga 
L 



tr{[G A lx (G K lx G K - G- lx G-) 



{G K lx G K - G_ lx G„) lx G R ]DG* } (A4) 



Here tr includes all summation in Dirac space and in- 
tegrations on space and time variables. We also de- 
note G- = G R — G A and index s refers to the relation 
Gg = 7oG if 7o- Also we use notation D = l/\z — y| and 
have dropped a term proportional to G R + Gf in the self- 
energy contribution which is zero by energy integration. 

The interaction part to the equilibrium noise S g = 
Se + S v can be calculated by applying the fluctuation 
dissipation theorem. Thus we have 



S g = ie ga 



(Pp d 2 q tr[I 2 (p)I 1 (q)+2J 2 {p)J 1 {q)} 
2tt (2tt) 2 |p-q| 



here the trace is taken over Dirac matrices. The functions 
Ji{p) and hip) include energy integration 



h{q) 

Hp) 

Ji(q) 
Hp) 



J ^fo(qo)G s -(qq ) (A5) 
- J f^(l - fZ( Po ))(G' lx G- lx G- + 



C. lx G. lx G') 



PPO 



(A6) 



-7o 



2tt 



(l-/ 2 (<Z ))(G_7xG_) g - go (A7) 



To 



/ ^ /o(Po)(G -1*G R + G A lx G _)(A8) 



here G' = ReG R . 

A direct energy integration yields 



and the vertex part acquires a form 

= ^^PYl I Tr["/oG a (yt 1 x.'t') lx T x G a (x't , zt 1 )io 
-(T x G <T (ztix.t)j x T x G <T (-x.tyt 1 ) + 

(A2) 



h(q) 
hip) 
Jlip) 

Mq) 



-ifo{q)i ■ q/<? 

16ttTA 2 (p)7 • p/p 2 

-8^707, A 2 (p) 



(A9) 
(A10) 
(All) 



2<f 



70(7 • q7*7 -<i + q ix) (A12) 



G a (zt 1 yLt)-i x T x G a (-x.tyt 1 )T x )] 



where symbol J denotes multiple integrations J = 
J dRd 2 rd 2 yd 2 zdh. 

The gauge transformation (fTTj) yields a V(t — t') de- 
pendence for all Green's functions. We simplify the noise 
by adding or subtracting expressions such as the energy 
integrals of the products of GF: G A G A or G R G R (which 
is zero). Thus for the zero frequency contributions to the 
noise power (at ui — > ) we get 



S v 



2e 2 iga 



tr[ 10 {G K lx G H + G A lx G K ) 
Div(G K lx G K -G„ lx G-)] 



(A3) 



where A(p) = <5/(7r(<5 2 + p 2 )) and q,p = |q|,|p|. Inte- 
grating over momenta in the formula for noise and col- 
lecting all contributions we note that logarithmic diver- 
gences are exactly compensated in the sum of two types of 
diagrams. Thus the Coulomb contribution to the equilib- 
rium (V = 0) noise becomes S g /A = 0.l9aTa m i n g. The 
prefactor 0.19 cames from numerical integration of the 
second term in S g and corresponds to the nondivergent 
part of vertex diagram (Fig. lb). 



Appendix B 

In this section we estimate two loop diagrams shown 
in Fig 2. The summations over the Keldysh indices leads 
actually, to the two sets of such diagrams: one is pro- 
portional to the sum of retarded and advanced Coulomb 
Green's functions, while the other is proportional to 
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Kcldysh Coulomb Green's function. As our principal ap- 
proximation we consider unscreened coulomb interaction 
to first order in g, while w-dependent parts of Green's 
functions include random phase approximation polariza- 
tion and, therefore, are a higher order in g. Thus, the set 
of diagrams that involves the Keldysh Coulomb Green's 
function is the second order in coupling constant, i.e. be- 
yond the 1st order considered here. As to the former 
diagrams, the relevant contribution may be expected in 
the first order in g. However, it simple show, and we 
proved this by direct calculations, that this set of dia- 
gram is proportional the small transmitted momentum 
q (in our case q — > 0) and does not contribute to the 
noise power. Indeed, in the limit of unscreened coulomb 
interaction (D R + D A )= l/\x — y\ the contribution to the 
noise which originates from this set of diagrams can be 
written as 



e 2 iga 
2L 



P{xtyt 1 ) = ^TrG (T (xty<i)7or x G (T (ytixt)7 x T x 



In equilibrium for Si oop we have 



Sioo P (q ->• 0,u) ~ iga min aP(qu)j^Q(qu) 



where we also introduced the Fourier transform for all 
functions. At zero temperature T = the current-density 
correlator P(quj) and the function Q(quj) can be easily 
calculated. In the limit of low frequency and uj > q we 
obtain 



dydzdt\P(xtyt\ 



1 



rQ(zt±X t ) 



P(quj) = 



Q{quj) 



V 1 - 1 2 /" 2 

-4iq x 



(Bl) 



-Po 



where Q(zt\x't') represents the current response to ex- 
ternal field and P(xtyt\) is the current-density correlator 
of nonintcracting electrons, 



Q(ztix't') 



TrG tT (zt 1 -x.'t')j x T x G a (-x't'zt 1 )j 



here q — \q\ and Po is the high momentum cutoff. Thus 
the shot noise contribution of this diagram for low fre- 
quency and uj > q is proportional the the small transmit- 
ted momentum q, i.e. vanishes as q — > 0. 
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